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■ ABSTRACT 

^ il It has been a long-standing question whether fossil groups are just sampling the tail 

of the distribution of ordinary groups, or whether they are a physically distinct class 
of objects, characterized by an unusual and special formation history. To study this 
question, we here investigate fossil groups identified in the hydrodynamical simulations 

■ of the GIMIC project, which consists of resimulations of five regions in the Millennium 

Simulation (MS) that are characterized by different large-scale densities, ranging from 
a deep void to a proto-cluster region. For comparison, we also consider semi-analytic 
models built on top of the MS, as well as a conditional luminosity function approach. 
We identify galaxies in the GIMIC simulations as groups of stars and use a spectral 
synthesis code to derive their optical properties. The X-ray luminosity of the groups 

' is estimated in terms of the thermal bremsstrahlung emission of the gas in the host 

\ halos, neglecting metallicity effects. We focus on comparing the properties of fossil 

groups in the theoretical models and observational results, highlighting the differences 
between them, and trying to identify possible dependencies on environment for which 
our approach is particularly well set-up. We find that the optical fossil fraction in 
all of our theoretical models declines with increasing halo mass, and there is no clear 
environmental dependence. Combining the optical and X-ray selection criteria for fossil 
groups, the halo mass dependence of the fossil groups seen in optical vanishes. Over 
the GIMIC halo mass range we resolve best, 9.0 x 10^^ ~ 4.0 x 10^^/i~^ Mq , the central 
galaxies in the fossil groups show similar properties as those in ordinary groups, in 
. terms of age, metallicity, color, concentration, and mass-to-light ratio. And finally, the 

satellite galaxy number distribution of fossil groups is consistent with that of non fossil 
groups. These results support an interpretation of fossil groups as transient phases in 
the evolution of ordinary galaxy groups rather than forming a physically distinct class 
of objects. 
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1 INTRODUCTION 

Galaxy formation is one of the most complex processes in 
the evolution of the Universe, and many aspects of it are 
still poorly understood. Based on the hierarchical formation 
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hypothesis, galaxy clusters and groups form through the as- 
sembly of many galaxies into a large common dark matter 
halo. Information about the merging history within a group 
is carried mainly by its primary central galaxy and the satel- 
lite galaxies. If evolved in isolation for a long time, the end- 
product of the merging in a group should be a system with a 
central early-type galaxy, and few left-over satellites. Groups 
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with a strongly dominating central early-type galaxy, and a 
bright extended X-ray halo with a cooling time scale longer 
than the group merging timescale, have been called "fossil 
groups" - because such an outcome could be most easily 
understood if these are groups with a very early assembly, 
leaving enough time for satellite galaxies to merge away 
with the centre. 

Aft er being or i ginall y discovered by IPonman et alj 

||1994 ), iJones et al.l (|2003l ) defined fossil groups (FGs) by 
requiring that these groups need to have a difference in ab- 
solute R-band magnitude of at least AAfi2 > 2 mag between 
the brightest and the second brightest galaxy located within 
half the projected virial radius, and that their X-ray lumi- 
nosity should exceed Lx.boi 2.13 x 10*^/i~^erg s~^. FGs 
have since been analyzed in numerous observational studies 
fe.p. Ijones et aL 2003; La Barbcra et al. 2009': I Santos et al.l 



I2OO7I : lMen~s de Oliveira et al., .200& ; .Aguerri et all l201l[) 



However, due to a rather low number of FGs identi- 
fied observationally and the lack of X-ray information 
in many cases, most works focused on the magnitude 
gap. This magnitude gap has been s tudied in theoretical 
work through large sky s urveys (e.g. van den Bosch et al] 
120071 : lYang et al.1 I2OO8I : iMilosavlievic et al.l 1200^) 7" pure 
dark matter simulations (e.g. iD'Onghia et al. l2007l : 
Ivon Benda-Bec kmann ct al. 20 0^ , semi-analyt ical models 



(e.g. 



, „ ISales et al.i2007l:lDi'az-Girnenez et al.ll2008l ') , or through 
gas simulations i Dariush et al.ll2007l 'l. However, these stud- 
ies do not agree well on the optical fossil fraction. 

A common theme driving many of these investiga- 
tions has been the desire to understand whether fos- 
sil groups are a clas s of objects that are special in 
their fo rmation history (iD'Onghia et aI]|2005l : [Dariush et al.l 



I2OO7I : I von Benda-Beckmann et al.l I2OO8I ). and whether 
they would show diff erent environmental dependences 
than ordinary groups l|yon Benda-Beckmann et al.l l2008l : 
iDi'az-Gimenez et al.l boilh . Some conflicting claims about 
this have been made, and it is not yet clear whether fossil 
groups can be easily accommodated in the leading ACDM 
model for cosmic structure formation. 

In this paper, we study the environment of fossil groups 
in the Galaxies Intergalactic Medium Interaction Calcula- 
tion (GIMIC) simulations, as well as in two semi-analytic 
models (SAM) of galaxy formation, and one conditional lu- 
minosity function (CLF) catalogue. Furthermore, we also 
investigate the properties of central galaxies in FGs and 
non-FGs through the GIMIC simulations. The GIMIC runs 
consist of re-simulations of five different density regions se- 
lected from the volu me of the Millennium Simulation (MS, 
ISpringel et al.|[2005l 'l. By identifying fossil groups in these 
regions, we can study in detail how the large-scale environ- 
ment impacts the properties of these systems, and for the 
first time, we can carry out such an analysis in full hydro- 
dynamic simulations that start directly from cosmological 
initial conditions. Furthermore, since the regions of GIMIC 
are drawn from the MS, we can read ily use the SAM cat- 
alogues of De Lucia fc Blaizoj (|2007l . hereafter D07), and 
I Bower et aL 1 20061 . hereafter BOG), to carry out a region- by- 
region comparison. In addition, the galaxy cata logues based 
on the conditional luminosity function model of lYangj l|201ll . 
hereafter CLF) have also been constructed for the MS halos; 
we can hence select the same five regions here as well. The 
comparison of these catalogues for the same regions yields 



insightful tests of the predictions of the theoretical models 
with respect to FG properties. We also compare the pre- 
dictions of these models with relevant observational data to 
assess how closely the models reproduce the observed FG 
properties. 

This paper is organized as follows. In Section [2j we in- 
troduce the GIMIC simulations, and describe our methods 
to analyze the runs in terms of a spectral synthesis code to 
compute galaxy properties and a simple approach to calcu- 
late the X-ray luminosity of host halos. In Section [S] the 
SAM and CLF catalogues we use are introduced. Section |3] 
is devoted to a description of the properties of the galaxy 
populations in the different model catalogues. A compari- 
son of the properties of FGs from these models is presented 
in SectionO along with their comparison with observational 
data. Finally, we discuss our results and summarize our main 
conclusions in Section 6. 



2 ANALYSIS OF THE GIMIC SIMULATIONS 

2.1 Simulation set and galaxy identification 

In the "Galaxies Intergalactic Medium Interaction Calcula- 
tion" (GIMIC) project, five different regions of comoving ra- 
dius ~20h'^ Mpc were selected from the volume of the Mil- 
lennium Simulation (|Springel et al.ll2005l ) and re-simulated 
with hydrodynamics included. The five regions have different 
mean over-densities that deviate by (— 2, — 1, +0, -1-1, -1-2) cr 
from the cosmic mean, where a is the rms mass fiuctuation 
on a scale of ~ 20 Mpc at z — 1.5. The simulations in- 
cluded physics modules for gas cooling and photoionization, 
quiescent star formation associated with supernovae feed- 
back, kinetic supernova feedback and chemodynamics, but 
no AGN feedback. The GIMIC simulations were run with 
GADGET-3, an updated ver sion of the GADGET-2 code (last 
described in 'Springell |2005| ) . with the physics implementa- 
tion corresponding to one of the models studied in the 'Over - 
whelmingly Large Simulations' project (|Schave et al.|[2O10l ). 
The GIMIC regions have been each evolved separately, and 
were run at several different resolutions. The mass of the gas 
particles in the GIMIC simulations was ~ IO^/i-^Mq, suffi- 
cient to make the galaxy catalog complete down to a r-band 
absolute magnitude limit of about —11.0. Throughout this 
paper, we adopt the same cosmological parameters as used in 
GIMIC and the parent MS project: flm = 0.25, Qa = 0.75, 
= 0.045, Us = 1, (78 = 0.9, Ho = 100 h km s'^ Mpc'^ 
with h = 0.73. A detailed description of further method- 
ological asp e cts of the GIMIC simulations can be found in 
ICrain et al.l (|2009l ). 

Halos in the simulation were identified by the stan- 
dard friends-of-friends (FOF) algorithm, applied to the high- 
resolution dark matter particles with a linking length of 
b — 0.2 in units of the mean inter-particle separation. Bary- 
onic particles were linked to groups by associating them 
with the group containing the nearest high-resolution dark 
matter particle. Substructures within FOF halos were then 
identified with the SUBFIND algorithm (|Springel et al.ll200ll : 
iDolag et al.ll200^ ). which decomposes each FOF group into 
a set of locally overdense and gravitationally bound struc- 
tures. For the larger FOF halos of group- or cluster-size, this 
procedure reliably identifies the different constituent galax- 
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ies that make up the groups. These appear as different sub- 
halos inside the FOF groups, and can be directly compared 
to observed galaxies. 

In optical observations, it has been found recently that 
not all the emission comes from stars that reside in the 
cluster's member galaxies. Instead, there is also a smoothly 
distributed stellar component which is typically peaked 
around the cluster's central galaxy, but extends to large 
radii. This diffuse stellar compo nent has been studied in a 
number of work s . For instance, iMurante et al.l (l2007h and 
iPuchwein etal\ l|2010l ) analysed high-resolution hydrody- 
namical simulations of galaxy clusters, and pointed out that 
this intra-cluster light (ICL) component accounts for up to 
~ 45% of all the light in simulated clusters and groups. 
This means we have to take this ICL component into con- 
sideration and can not simply add it to the luminosity of 
the central group galaxy, w h ich wo uld artificially boosts its 
luminosity. Puchwein et al.l (|201(]|) discussed four different 
methods to identify the ICL component, all of them showing 
quite similar results. So we here simply follow their simplest 
method, limiting the aperture, to exclude the ICL of central 
galaxies. The cut radius is given by 



rcut = 27.3 h ^kpc x 



crit200 



(1) 



where Mcrit200 is the mass within the radius enclosing a 
mean overdensity 200 times the critical density of the Uni- 
verse. In Figure[T] we show the measured luminosity function 
of galaxies in GIMIC, with (red solid line) and excluding 
(red dotted line) the ICL stars. The two luminosity func- 
tions are not significantly different. Unless stated otherwise, 
we will always use central galaxies without ICL component 
throughout this paper. 



2.3 X-ray luminosity of group halos 

Because one of the primary criteria often used to distinguish 
fossil groups from ordinary groups is the X-ray luminosity, 
we measure this quantity directly in the hydrodynamic sim- 
ulations for all of our identified groups. We adopt a simple 
approach to calculate the X-ray luminosities in the GIMIC 
simulations, where for simplicity we neglect metal lines and 
only consider the continuous X-ray emission that is due to 
thermal bremsstrahlung. Neglecting the contribution from 
metal lines can provide an underestimate of X-ray luminosi- 
ties, especially at relatively low X-ray temperatures, T < 3 
keV. This might undersestimate the bumber of X-ray se- 
lected fossil groups in our sample, but does not affect their 
fossil optical properties. We model the hot gas inside subha- 
los as only containing a primordial mixture of hydrogen and 
helium, assuming that it is completely ionized and optically 
thin. 

The bolometric X-ray luminosity Lx of a group halo 
can be obtained by integrating the thermal bremsstrahlung 
emissivity over the halo's volume. To calculate this integral, 
we can cast the volume integral into a discrete sum over the 
SPH particles. This gives 



dE 

dF = '^ 
0.864 ft: 



0.72gs T2 p2 



dV 



V 
N 



(3) 



Pi mi, 



where qb is the bolometric Gaunt factor, which we set to 
1.2, nip the proton mass, k = (^^)^^;^^: with me for 
the electron mass, e for the electron charge, and T for the 
local plasma temperature. In the above equation, the volume 
element associated with each of the particles is AVi = pi/rrii, 
where pi and irn are the SPH density and particle mass for 
the i-th particle, respectively. 



2.2 Optical galaxy properties 



We treat each stellar particle in the simulations as a simple 
stellar populati on with the initial mass function (IMF) of 
IChabrieJ (|2003l '). with a mass range of 0.1 — 100 M©, taking 
into account the age and metallicity of the star particle. 
This IMF matches the one that has also been used when 
running the GIMIC simulations. The spectrum of the stellar 
p opulation is produced by i nterpolating the SSP templates 
of iBruzual fc CharlotI l|2003l V 

The spectrum of each simulated galaxy is then simply 
obtained by adding up all the stellar particles' spectra within 
the same subhalo: 



itot ( A) = ^ rrii Lspp{ti, Zi, X), 



(2) 



where rm, ti, and Zi are the mass, age, and metallicity of 
stellar particle i, respectively. A'" is the total number of stel- 
lar particles within the subhalo, and Lspp is the intrinsic 
SED per unit mass interpolated from the SSP templates. 
Having obtained the spectra, we then apply the SDSS fil- 
ters in the AB system to get the galaxies' u, g, r, i, and 
z-band magnitudes. We do not account for dust extinction 
in our procedure. 



3 SEMI-ANALYTIC AND CONDITIONAL 
LUMINOSITY FUNCTION MODELS 

3.1 Semi-analytic galaxy catalogues 

The Millennium Simulation has been used extensively by 
the Munich and Durham groups to construct so-called semi- 
analytic galaxy formation models in which the evolution of 
galaxies is followed based on dark matter merger histories 
measured from the simulation. The baryonic physics is de- 
scribed at a coarse level with a simplified set of differential 
equations, allowing it to be integrated forward in time com- 
paratively quickly. The two semi-analytic galaxy formation 
codes differ in their detailed assumption about the param- 
eterization of the relevant physics, making it interesting to 
also compare them with each other, prompting us to con- 
sider both models. 

Full de tails of the Munich semi-anal y tic model can 
be found in Kauffmann fc Haehnelt (^2000^■ Springel et al.l 
ll200lD. iDe Lucia et al.l l|2004h . and (Pe Lucia et al.l (|2006l ). 
Croton et al.l ||2006|) . We here analyze the model by 



De Lucia fc Bla izoT ("20 071, D 7), w hich is itself a modi- 



fied version of Croton et al.l (|2006l ). and includes a pre- 
scription for the growth and activity of central black holes 
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Figure 1. The galaxy luminosity functions of the GIMIC simulations (solid red line), GIMIC simulations without IGL (dotted red 
lines), the GLF catalogues (solid magenta) of I YantJ ll201lh . and the two SAMs from B06 (dashed blue line) and DOT (solid magenta 
line) in the five different regions we analyzed, labeled as +2a to —2a fro m most overden s e to m ost underdense. In all panels, the black 
solid line shows the Schechter form luminosity fitting function obtained bv lBlanton et al ] l|200ll ') from the SDSS observational data. The 
lower-right panel shows the averaged luminosity function of the five regions. 



and their effect on suppressing the cooling and star forma- 
tion in massive halos. D07 used the initial mass function of 
IChabrieil l|2003l ). Their publicly available catalogues include 
26,787,155 galaxies at 2 = in the whole Millennium vol- 
ume. The halo mass in this catalogue is Mcrit200, which is 
defined as mass within the radius enclosing a mean overden- 
sity 200 times the critical density of the Universe. The dust 
extinction is included in SDSS observer frame magnitudes 
for this galaxy catalogue. 



The Durham semi-analytic galaxy formation model 
is implemented in the GALFORM code (ICole et al" 



iBenson et all l2003l : iBaugh et all l2005l : iBower et all 



The semi-analytic galaxy catalogue we analyze here is de- 
scribed in detail bv lBower et al. 1 (120061 . B06), and contains 
several enhancements compared to older GALFORM mod- 
els, such as the formation and growth of black holes, disc in- 
stabilities, improved cooling calculation and AGN feedback. 
The B06 model adopts a Kennicutt IMF with no correction 
for brown dwarf stars. The total number of galaxies number 
in the final z = catalogue is 24,537,589. In this catalogue, 
dust extinction is also included, but the magnitudes are in 



the Vega system, which we change to the AB system to 
match the other catalogues used in this paper. 



3.2 Conditional luminosity function catalogues 

The conditional luminosity function (CLF) technique is a 
statistical procedure to populate dark matter halos with 
galaxies such that the observed galaxy luminosity function 
(or stellar mass function) and the clustering properties of 
galaxies can be accura tely reprodu c ed. Th e method has been 
described in detail by lYang et al.l l|2003l ), and has bee n re- 
fined over the years in numerous studies (lYanell201ll . and 

references therein). 

In a recent work, lYangj i 201lD were able to take the 
full information of subhalo merger histories into account 
when constructing the CLF as a function of redshift. Here 
we make use of the model parameters that are constrained 
using the SDSS luminosity (and stellar mass) function and 
the projected two point correlation functions. MS dark mat- 
ter halos are populated with galaxies in two catalogues: one 
constructed with galaxy luminosities, the other with galaxy 
stellar masses. 
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In the luminosity catalogue, each galaxy is assigned an 
r-band absolute magnitude °'^ Mr — 5 log h, which i s K + E 
corrected to redshift 2 = 0.1 l|Blanton et al.ll2003h . While 
in the stellar mass catalogue, the stellar mass is assigned so 
that we can repro duce the stellar mass function obtained by 
lYang et alj l|2009l ). 

Note that as tested in the CLF models and 

thus the constructed galaxy catalogues can not only repro- 
duce the observed SDSS luminosity (stellar mass) functions, 
but also halo occupation numbers obtained from the SDSS 
group catalogues. These catalogues contain roughly a total 
of lO'^ galaxies, slightly less than the two SAM models be- 
cause of the brighter luminosity (stellar mass) cut. 



3.3 Matching of the GIMIC regions and 
comparing the four models 

The spatial location and radius of the five different density 
regions re-s i mulate d in GIMIC can be found in Table 1 of 
ICrain et al.1 (120091) . It is hence straightforward to cull the 
semi-analytic models and the CLF catalogues to the same 
regions. Also, all the models have a parameter that identi- 
fies whether a galaxy is the central galaxy in a group, or 
whether it is a satellite galaxy within a group. We test only 
the central galaxies for membership in one of the regions, 
and if this is the case, always use all the galaxies in the 
same group, too. 

The number of galaxies with "'^Mr — 5 log h < —IQ and 
the number of FOF groups (in brackets) in each region is 
listed in Table [1] The four models do not agree with each 
other well in the number of galaxies and the number of FOF 
groups in all five regions. The reason is that the four models 
have different predictions for galaxies even in the same dark 
matter halos. With the magnitude limit we applied, part of 
low mass galaxies and/or groups are not included. It is inter- 
esting to see that the GIMIC model produces fewer galaxies 
compared to other models in all regions for this magnitude 
cut, especially in the +2a region. From Figure [T] the drop 
of $ in GIMIC is mainly caused by the strong wind model. 

Note that since the absolute magnitudes provided in 
SAMs and GIMIC are scaled to redshift z = 0, in or- 
der to properly compare with the SDSS observation, we 
have converted them into the ones corresponding to redshift 
z = 0.1 using the av e rage K + E correc t ions p rovided by 
iBlanton et aLl l|2003D : iBlanton fc RoweisI l|200'<i) . Through- 
out the paper we use the same "'^Mr — 5 log h notation. The 
four models also adopt different definitions of halo masses. 
For proper comparisons, we convert the halo masses in the 
SAMs and GIMIC to the ones in the CLF model with a halo 
abundance matching method. The final halo masses thus 
correspond t o the FOF halos w hose mass function is well 
described bv lSheth et al.l (|200ll ). We apply this normalized 
halo mass of each model in the following sections. 



4 BASIC SAMPLE PROPERTIES 

Before we analyze the specific properties of fossil groups, we 
study here some of the basic properties of our different model 
galaxy catalogues, and how they compare to observational 
data. We limit ourselves to a discussion of the luminosity 



function, the properties of central galaxies in groups, and 
the halo occupation distribution function. 

4.1 Luminosity function 

The luminosity function is one of the most fundamental 
characteristics of a galaxy population. In Figure [TJ we show 
the GIMIC luminosity function of the five different regions 
in separ ate panels, comparin g with the observational re- 
sults of IBlanton et al.l (|2001l ) and with our semi-analytic 
and CLF models. A Schechter function fit provides a good 
description of the dust-corrected SDSS data over a range 
—23 < "'^Mr — 51og/i < —16, with best-fit parameters 
0* = 1.46 X 10"^/l^Mpc"^ Af* = -20.83, and a = -1.20. 
The luminosity functions in the different density regions are 
clearly different in their normalization. The higher density 
regions have a higher number density of galaxies, as ex- 
pected. However, there is no obvious difference in the shape 
or characteristic magnitude. In the lower-right panel of Fig- 
ure[T] the total luminosity functions from the CLF and SAM 
models are also shown; they exhibit a better mat ch than 
GIMIC to the Schechter fit of IBlanton et al.1 (120011 ) for the 
observational data. We note that the CLF and SAM lumi- 
nosity functions have a higher normalization than the obser- 
vations. This is mainly due to the fact that we have added 
all the galaxies in the five regions with equal weight, for sim- 
plicity, while in reality there are fewer high density regions 
than low density regions within a representative volume of 
the Universe. 

It is obvious from Figure [1] that there is a substantial 
disagreement between the shape of the LF of the GIMIC 
simulations and the other theoretical models. In fact, the 
concave shape and bump- like feature in the GIMIC luminos- 
ity function appears in all different density regions, and is 
caused by the specific feedback mo del that was i ntrod uced in 
the hydrodynamic simulations. In lCrain et al.l ((200^ it was 
shown that the galaxy stellar mass fun ction has a simila r 
trend (see the right panel of figure 3 in ICrain et al.ll2009l ). 
They also investigated the origin of the disagreement in some 
detail. At the lower mass end (M < 10^h~^ Mq), the excess 
number of galaxies reflects a reduction in the efficiency of 
SNe feedback in poorly resolved galaxies, while at interme- 
diate masses (10^ ~ W^'^h~^ M©), the 'dip' is caused by an 
efficiency peak in the feedback wind model, that is even too 
effective in quenching star formation for galaxies in this mass 
range. Finally, at the high mass end (M. > 10"/i"^ Mq), 
the overprediction of massive galaxies is mainly caused by 
gas overcooling in the brightest cluster galaxies, which is 
not efficiently counteracted by SN -driven feedback (see also 
ISaro et al.ll2006l : iNuza et all I2OI0I ') . Resolving this problem 
may require the introduction of AGN feedback. 

Taken at face value, the poorer agreement of the hy- 
drodynamic GIMIC simulations with the observational data 
may seem to suggest that these calculations are inferior to 
the other theoretical models. It needs to be taken into ac- 
count, however, that the hydrodynamic simulations do not 
involve much tuning, intead they yield a relatively precise 
and physically self-consistent calculation of the outcome of 
cosmic evolution for a specific set of assumptions about star 
formation and its regulation through feedback. The simu- 
lations show that it is actually very hard to come up with 
a successful parameterization of the feedback processes, ar- 
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Figure 2. Upper panels: Median r-band luminosity of the brightest cluster galaxy (BCG) as function of halo mass of four theoretical 
models, as labeled. Lower panels: Median stellar mass of the most massive galaxy (MMG) as function of halo mass of the models. The 
solid lines in each panel give the results for the first ranked galaxies in the groups. The different colors and line styles encode the different 
density reg i ons, a ccording to the key shown in the lower-right panel. For comparison, the solid black lin e shows the fitting formula of 
lYang et al] l l200d) derived for the CLF model. The shaded area and error bars come from I Yang et all ||2008| ). and show the 68% confidence 
regions of halo mass estimated through the group luminosity Mj^ and group stellar mass Mg, respectively. 
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Figure 3. Upper panels: The luminosity ratio between first and second ranked galaxies, as a function of halo mass for our four theoretical 
models, as labeled. Lower panels: The stellar mass ratio between first and second ranked galaxies as a function of halo mass for the 
same models. The colored lines encode the different density regions, using the same key as in Figure [2] The cr oss and square sy mbols 
represent an analysis of the observa t ional data, only differing in the treatment of the second brightest galaxy l|Yang et al.ll200a '). The 
shaded areas come from pYang et af] l|2008h . which indicate the corresponding 68% confidence intervals for crosses (red shaded region) 
and squares (blue shaded region). 
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Model 


+2a 


+ la 


Oct 


-Icr 


-2a 


GIMIC 
D07 
B06 
CLF 


4,732 (2,583) 
14,516 (7,386) 
10,217 (5,298) 
13,553 (5,256) 


1,418 (809) 
3,576 (2,226) 
2,431 (1,512) 
2,188 (1,610) 


1,046 (651) 
2,394 (1,500) 

1,547 (961) 
1,437 (1,060) 


898 (678) 
1,957 (1,525) 
1,212 (958) 
1,303 (1,115) 


591 (471) 
1,261 (1,022) 
736 (620) 
822 (732) 



Table 1. The raw number of galaxies with ''■^M,- — 5 log h < —16 in the five regions in each of the different theoretical galaxy catalogues. 
For reference, the second numbers given in brackets list the number of FOF-groups to which the galaxies belong. 



guably harder than it may seem based on the success of 
the semi-analytic models. This clearly indicates that the ef- 
fect of relevant feedback processes is still poorly understood 
in direct hydrodynamical simulations. However, one should 
also remember that the success of SAMs in reproducing lu- 
minosity function data is also the result of a tuning of the 
parameters entering the models. 

4.2 The properties of central galaxies 

To investigate the suitability of GIMIC, SAM and CLF 
galaxy populations for fossil groups, we next consider prop- 
erties of the brightest cluster galaxy (BCG) and the second 
brightest galaxy. Yang et al.l (|2008 ). based on the galaxy 
groups extracted from SDSS DR4, produced two fitting for- 
mulae for the relation b etween the lumin osity of BCGs and 
the halo mass (Eq. 6 in I Yang et al.|[2008l ). as well as for the 
relation of stellar mass of the most massive gala xy (MMG) 
and the halo mass (Eq. 7 in I Yang et al]|2008l ). The halo 
masses in their group catalogues, and thus in the two fitting 
formulae, are estimated though the ranking of the character- 
istic group luminosity Ml and the ranking of th e character- 
istic group stellar mass Ms (see more details in lYang et al.l 
The two sets of halo masses agree reasonably well 
with each other. Note however that the halo masses in their 
models are obtained for the WMAP3 cosmology. For con- 
sistency, we here converted the halo masses in their models 
to the ones corresponding to the MS using the abundance 
matching method. We use these converted fitting formulae 
as a reference standard. The different halo masses of each 
model are normalized to the CLF model, as described in 
Section |3] The selection of th e central galaxy da ta from our 
simulation is the same as in lYang et al.l (|2008l ): either the 
brightest cluster galaxy (BCG) or the most massive cluster 
galaxy (MCG) is taken as the central galaxy. For the "iso- 
lated galaxies", we apply a cut with the magnitude limit 
of "■'^Mr — 51og/i < —16. To change to a common mag- 
nitude system, we adopted for all models an r-band mag- 
nitude equal to r = 4.7 6 for the Sun in the AB system 
(jBlanton fc RoweidlioOTl ). 

In the upper panels of Figure [51 we compare model pre- 
dictions and observational data for the relation between the 
r-band luminosity of the BCG and the halo mass. In a sim- 
ilar way, in the lower panels we show the relation between 
stellar mass of the MMG and halo mass. Different panels are 
for the different model galaxy catalogues. In each panel, re- 
sults are shown for the regions at different overdensities. The 
SAM models and the CLF model both follow the observed 
relations reasonably well, while the GIMIC simulations show 
substantial disagreements. The tendency of GIMIC to over- 
predict the mass of the MMG and the luminosity of the BCG 



begins at halo masses of ~ 10^^ Mq. This overprediction 
occurs because GIMIC forms too many star particles in the 
central galaxies at late times. 

We now explore the gap between the first and second 
ranked galaxies. In Figure [3] we show the ratio between the 
brightest (massive) and the second brightest galaxy versus 
the mass of the host halo, which is directly related to the fos- 
sil characteristics of groups. The cross and diam ond symbols 
inclu ded for comparison come from the SDSS (lYang et al.l 
I2OO8I ) , and correspond to two different treatments adopted in 
case the second galaxy in a group is not observed, because it 
is fainter than the SDSS magnitude limit of about r ~ 17.7. 
In the first case, the second brightest (massive) galaxy is 
treated as having zero luminosity (mass), and is shown with 
crosses. In the second case, the second brightest (massive) 
galaxy is treated as the magnitude limit of the SDSS sur- 
vey (diamonds). The two SAM models and the CLF model 
all show a better fit to the observational data than GIMIC 
at halo masses greater than 10^^/i~^ M0. In lower density 
regions, the ratio in the CLF model and in the two SAMs 
also does not show large deviations from the higher density 
regions, unlike GIMIC; but here no tight observational con- 
straints are available. As the observation prefered, the more 
massive haloes tend to have smaller luminosity gaps. This 
is supported b y our cata logues, and also consistent with the 
findings of ,Skibba et all (|2007..,2011i ). 



4.3 Halo occupation distribution 

In Figure 131 we plot the satellite numbers versus halo mass, 
with a specified magnitude limit, "'^Mr- — 5 log h < —19.0 , 
for aU models. The fitting formula of lYang et all (|2008l ) 
that we include for comparison takes a power law form, 
{Ns} = (Mh/Msfi)'' , where (Na) is the mean number of 
satellite galaxies, Mh is the halo mass, and the two param- 
eters A/3,0 and 7 have the best fitting values of 10^'^ '^'^ Mq 
and 1.06, respectively. After accounting for the cosmology 
transformation from WMAP3 to Millennium for the fitting 
formula, and a renormalization of the halo mass to CLF 
using halo abundance matching method for all the models, 
the satellite numbers of all the models are roughly consistent 
with the observed power-law fitting relation from observa- 
tional data, except at the low halo mass end. Quite inter- 
estingly, based on the results shown in Figure |4l we do not 
see any evidence for environmental dependence of the halo 
occupation distribution in any of the considered models. 
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Figure 5. The optical fossil fraction as a function of halo mass for all four theoretical models we analyzed here. Each panel shows the 
results for one of the catalogues, as labeled. The lines with different colors and styles correspond to the different density regions, as in 
Fig.H 
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Figure 6. Left panel: The optical fossil fraction of GIMIC as a function of halos mass. Middle panel: The fraction of GIMIC groups 
with X-ray luminosity Lx larger than 2.13 X 10*^ h~^ergs~^. Right panel: The overall fossil fraction when both selection criteria are 
combined. In each panel, we show results for different overdensities using different colors and styles, based on the same key as in Fig. [2] 



5 RESULTS FOR FOSSIL GROUPS 

5.1 Abundance of fossil groups in different 
environments 

One of the most important parameters for fossil groups is 
their abundance. Previous analyses have studied, for ex- 
ample, the 2-deRree Field Galaxy Redshift S urvey (e.g. 
van den Bosch et all 120071 ') . the SDSS -DR4 jYang et al.l 



20081 ). semi-analvtic galaxy catalog ues (ISales et al 



or cosmological N-body simulations (|D'Onghia et al.| 



2007 



2007 



Since X-ray data has often been unavailable, many works 
only considered the magnitude gap between the first and 
second brightest galaxies, disregarding the requirement 
for a group to be also X-ray luminous. We note that 
the predicted optical fossil fractions from these works do 
not agree particularly well with each other, demonstrat- 
ing that the uncertainty in this quanti ty is still high 
(see Table [21). Interestin gly, seve ral studies ( [ P'Onghia et al.l 
2005; Milos avlieyic et al. 2006 



i 



van deiL Bosch et al 
2010t l 



20071 : 



Yang et all l2008l : iDariush et all |2007| . I2010f l have consis- 
tently claimed a trend of an increasing optical fossil fraction 
with decreasing halo mass. 

We start by defining fossil groups based only on the 
magnitude gap in the r-band, requiring it to be larger than 2 
for an optical fossil group, which is the commonly employed 



definition. In Figure (5] we plot the resulting average optical 
fossil fraction as a function of halo mass. A mass trend is 
clearly present in all models. It is not surprising to see such 
a declining trend in light of what we showed in Figure O 
The strong disagreement in GIMIC (magenta line) is pri- 
marily related to the sparseness of the data for the highest 
halo mass region, and the associated statistical uncertainty. 
However, the optical fossil fractions in each model are not 
fully consistent in detail. The error bars in Figure [5] reflect 
the la scatter obtained from 200 bootstrap resamplings of 
different groups. All models predict a very high optical fossil 
fraction at the lower halo mass end. This means that large 
magnitude gaps become much more common in low mass 
halos. Taken into account the big error bars, we do not find 
any clear evidence that environment has an effect on the 
optical fossil fraction in any of the models. 

As we discussed earlier, the hot gas profiles in the 
GIMIC halos can be used to measure the X-ray luminos- 
ity Lx of halos, providing us with a way to arrive at higher 
specificity in the selection of fossil groups. We now define a 
group to be fossil only if its X-ray luminosity Lx is larger 
than 2.13 x 10*^ /i~^ergs~^, in addition to showing the 
magnitude gap. This X-ray luminosity limit first adopted 
by [jpnos ot alj (|2003D is meant to exclude normal galax- 
ies which are not elliptical galaxies at the centre of groups. 
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Mass range {h~^ Mq) 


Optical fossil fraction 


Fossil fraction 


Reference 


1013 - 1014 


6.5% 


- 


van den Bosch et al. (2007') 


10^2 - 10i3 


13.4% 


- 


van den Bosch et al. (2007) 


~ 1013-5 


11% - 20% 


- 


Yang et al. f2008~) 


~ 1013 


18% - 60% 


- 


Yang et al. f2008~) 


1 \y 1 nl3 


18 /c 




U Ungnia et ai. izUU/ ) 


^ 1013 


8-10 




Sales et al. C2007~) 


~ 10l3 - IQlS 


13.3% 




Dariush et al. (2007) 






7.2% 


Dariush et al. (2007) 


^ 1013 _ 10" 


5 - 40% 




Milosavlievic et al. (2006) 






8 - 20% 


Jones et al. ("2003") 


~ 1014 


33% 




D'Onghia et al. (2005): Sommer-Larsen (2006) 



Table 2. Compilation of fossil fractions reported in the literature. In the first column we list the halo mass range for which the estimate 
applies, and where possible, we distinguish between optical and X-ray based fossil fractions. The last column gives the corresponding 
references, which here include both theoretical and observational works. 
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Figure 4. Mean satellite number as a function of halo mass. The 
colored lines encode the different density regions in the same way 
as in Fig. [2] for an absolute magnitude limit '^ iM^ — 5 log h < 
— 19.0 . The solid black line is the fitting result by lYang et al.l 
for the same magnitude limit. 



and inevitably imposes a selection bias on the data. In Fig- 
ure [6l we show the fossil fraction in GIMIC as a function 
of halo mass based on this definition. In the middle panel 
of figure [6l we show the fraction of groups that have X-ray 
luminosity Lx above the threshold 2.13 x 10*^ /i~^ergs~i. It 
is not surprising to see that the fractions show an increasing 
trend with halo mass. Although there are more optical fossil 
groups in low mass halos, the fossil fraction is very small 
(less than 15%) due to the limited X-ray luminosity. In gen- 
eral, because of the additional X-ray selection, the halo mass 
dependence of the fossil groups seen in the optical band van- 
ishes. 

Figure [S] also reveals that X-ray bright groups in low 
density regions are very rare, which appears consistent with 




40 42 44 

L, [h"' erg s"'] 

Figure 7. The group abundance as a function of X-ray luminos- 
ity. The quantity n shown on the y-axis is the group number in 
each Lx bin, while N is the total group number. Differently col- 
ored and styles lines represent the different density regions, as in 
Fig.H 



the dramatic drop found for low mass halos bv lDariush et al.] 
(2007, sec their Figure 5). To investigate this reduction of 
the X-ray fossil fraction in two lower density regions, we 
further plot the abundance of the groups in Lx-bins for the 
halo mass range IQi^ - W^''h~^ Mq in Figure|7| Most of the 
X-ray luminous groups in the two lower density regions are 
located around Lx ~ 10^* h~^ergs~^ . There are however 
very rare galaxy groups with Lx > 2.13 x 10**^ h~^evgs~^ 
in the two lowest density regions. The abundance of systems 
with Lx > 2.13 X 10*^ /i~^ergs~i within the average den- 
sity and high density regions are much higher than in the low 
density regions. The main reason for this behaviour is that 
structures in low-density regions are generally characterized 
by shallower potential wells. Therefore, the gas reaches rela- 
tively lower densities in the central part of these structures, 
thus explaining the correspondingly lower X-ray luminosity. 



5.2 The properties of fossil groups 

We now focus on the halo mass range 9.0 x lOi^ ~ 4.0 x 
lO^^/i"^ Mq, which is likely the most reliable in GIMIC. In 
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Figure 8. Basic central galaxy properties of fossil and non-fossil groups. The individual panels give the relations between six quantities: 
the X-ray luminosity Lx, the r band magnitude Mr, the mean age in Gyrs, the mean metallicity Z, the {g — i) color, and the galaxy 
concentration C. Crosses represent non-fossil groups, while solid circles are for fossil groups. The different colors refer to the different 
density regions as in our previous figures. 



Figure [S] we summarize our measurements for some of the 
basic properties of the central galaxies of fossil groups and 
non-fossil groups, such as age, metallicity, colors or concen- 
tration, and we discuss them in the following. 

IJones et al.l (|2003l ) found that fossil groups are more 
X-ray luminous for a given optical luminosity than non- 
fossil groups. Using Chandra X-ray observations of five fossil 
groups, suppleme nted by additional systems from the liter- 
ature, iKhosrosha hi et al. confirmed their result of a 
Lx-Le, relation. However, a recent work bv lVoevodkin et al.l 



(|2010t ) reached a different c o nclusi on. T hey claim that the 
compa risons in Ijones et al.l (|2003l ) and IKhosroshahi et al] 
were affected by systematic errors due to a non- 
uniform selection, because the data for fossil groups and 
other systems were obtained separat ely. Using instead the 
data from a uniform sample, iVoevodkin et al.l (j2010) con- 
cluded that fossil groups are not systematically brighter 
in X-rays than other groups at the same optical luminos- 
ity. Our results in Fig. [5] for the M r -Lx relation support 
the conclusion by IVoevodkin et al] l|201Cf ). In the region 
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Lx > 2.13 X W^^ /i~^ergs~^, we do not see a systematic 
increase of the X-ray luminosity for fossil groups at a given 
optical magnitude. 

To quantify the age of our simulated galaxies, we di- 
rectly use the mean age of the stellar particles. Our results 
for this analysis are shown in the third row of Figure [S] Fos- 
sil groups in these plots typically have a mean age around 
9 Gyr, but this age is not significantly enlarged compared to 
the age of non-fossil groups. This suggests that there may 
be little justification in the end to call the magnitude-gap 
selected sample of groups "fossil" - they do not appear to 
have different ages than normal groups. 

Similarly, the mean galaxy metallicities for fossil and 
non-fossil groups reported in the fourth row of Figure[H]show 
no apparent differences. We have also used our spectral syn- 
thesis code to calculate colors for the simulated galaxies. 
The g~i color is shown in the fifth row of Figure [S] We here 
find values that are systematically lower than observational 
results. That is because we have not introduced any dust 
attenuation model in our synthesis code, and also because 
the simulation tends to have too much recent star formation 
in central galaxies. This is also the reason why the colors of 
our simulated fossil groups spread over a large range from 
0.5 to 1.0 in the plots. Many of them are actually too blue 
due to recent star formation to represent bona-fide early 
type galaxies, even though their large stellar mass otherwise 
suggests that they are. 

Finally, the last row in Figure[5]shows the concentration 
C of the central galaxies. We here used the maximum circu- 
lar velocity Knax of the central galaxy's halo and the radius 
''max at whi ch this velocity is att ained to estimate halo con- 
centrations jSpringel et al.llioosi '). Note that this procedure 
assumes that the dark matter halo is not substantially modi- 
fied by the settling of the baryons in the center, and that the 
rotation curve is dominated by dark matter at its maximum; 
both of these assumptions may actually break down if the 
central galaxies are very luminous and compact. Regardless 
of this complication, we may still go ahead and compare the 
concentration estimates obtained in this way among the dif- 
ferent regions, as done in Fig. [8] We see that although there 
is one fossil galaxy in the +2a region that shows an un- 
usually high concentration value, most fossil groups are well 
intermixed with the distribution of normal groups, providing 
no evidence for a significant difference in the concentrations. 

We have also investigated the mass-to-light ratios of 
the central gala xies in the selected GIM I C data. Previ- 
ous works (e.g. Khosrosh ahi et al.l l2007l : IVikhlinin et al.l 
I1999I : lYoshioka et a l. 2004 1 have found conflicting results 
for the mass-to - light ratio of fossil groups. Using M/Lr 
IVikhlinin et al.l (| 19991 ) reported a mass-t o-light ratio ~ 
250 — 450Mq/Lq, and using Mcrit2oo/ Lb lYoshioka et al.l 
(|2004l ) found ~ 100 - IOOOMq/Lo. Both of these stud- 
ies claimed a very high mass-to-light ratio for the X- 
ray overluminous elliptic al galaxies, while using Mvir/is, 
iKhosroshahi et al] (|2007l ') reported that the mass-to- light ra- 
tio of fossils tends to be at the upper envelope of the values 
seen in normal groups and clusters. In Figure [9l we show our 
measurements for the GIMIC simulation. We use Afcrit200 as 
the group mass to calculate the mass-to-light ratio. Here we 
also find that the mass-to-light ratio for fossil groups is con- 
sistent with that found for non-fossil groups. 

Based on our results, there is hence little tangible ev- 




13.2 



Figure 9. Mass-to-light ratio in the GIMIC simulations as a 
function of galaxy stellar mass. The crosses are central galaxy's 
mass-to-light ratio from non-fossil groups, while the solid circles 
show fossil groups. The different colors encode the different den- 
sity regions, as in the previous figures. 



idence for differe nces in the prope r ties o f fossil and non- 
fossil groups. In iLa Barbera et al] (|2009l ). they compared 
both structural and stellar population properties of fossil 
and non-fossil group galaxies, which is also shown consistent 
results. This calls into question whether the fossil groups 
are indeed a special class of objects that is characterized 
by an unusual formation history. Instead, our results appear 
more consistent with an interpretation of fossil groups as a 
common phase in galaxy evolution during which they lie in 
the tail of the distribution of a certain group properties, in 
this case the magnitude difference between first and second 
ra nked galaxies. This re sult is in line with findings obtained 
bv lDariush et al.l (|2010l ). who also used data from the Mil- 
lennium simulation. These authors found that, no matter at 
what redshift the fossil groups are selected, after ~ 4Gyr 
more than ~ 90 percent of them change their status and 
become non- fossil according to the magnitude gap criterion. 
This behavior, obtained from the galaxy formation histories, 
directly supports our conclusion here. As a further check, we 
show in Figure [TU] the number of satellites for fossil and non 
fossil groups. Here we use all the satellite galaxies in the 
GIMIC simulations. The satellite number for fossil groups 
is not different from that obtained considering non fossil 
groups, which confirms that the fossil groups do not have 
unusual formation histories. Nevertheless, we stress that al- 
though the fossil groups defined according to the magnitude 
gap are statistically similar to non fossil groups, there may 
be rare situations in which the fossil groups indeed have 
unusual formation histories (the few outliers visible in Fig- 
ure [lOl). 
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Figure 10. The number distribution of all satellite galaxies in 
the GIMIC simulations. The different lines are the mean satellite 
number for non fossil groups, while the diamonds show results for 
fossil groups. The different colors correspond to different density 
regions, as in the previous figures. The error bars indicate the 
r.m.s. scatter within each mass bin. 

6 DISCUSSION AND CONCLUSIONS 

In this study, we have analyzed fossil groups in the hydrody- 
namic simulation data of GIMIC combined with CLF cata- 
logues and two semi-analytic galaxy formation models. Our 
primary goal has been to shed light on the question of how 
well models of galaxy formation reproduce the observational 
properties of fossil groups, and whether there is any evidence 
that this class of objects is characterized by a special forma- 
tion history. To this end we have first investigated the lumi- 
nosity functions of our theoretical models, and in particular 
the abundance of brightest and most massive cluster galax- 
ies. We then considered the ratio of luminosity or stellar 
mass between the first and second ranked cluster galaxies, 
and examined the satellite abundance as a function of host 
halo mass. In all of our measurements, we have also checked 
whether there is any environmental difference between the 
five regions simulated by GIMIC, which have substantially 
different mean densities. 

For defining "fossil groups", we have mostly employed 
the luminosity gap between the two first ranked galaxies 
in a group, similar to what has been done in many pre- 
vious works. However, by also measuring the X-ray lumi- 
nosity of halos in GIMIC, we could more sharply define 
fossil groups as those which also have high X-ray lumi- 
nosity, refiecting more closely the originally used observa- 
tional definition. Since X-ray emission is unavailable in the 
CLF and the SAMs, we have however applied this selection 
only in parts of our analysis. Finally, we have investigated 
some basic properties of fossil and non-fossil groups in or- 
der to see in which properties they differ most significantly. 
For this, we restricted our halo sample to the mass range 
9.0 X 10^^ ~ 4.0 X 10"/i"^ M0, which should be the most 
reliable in our hydrodynamical simulations. 
Our primary conclusions can be summarized as follows. 

• The optical fossil fraction in all of our theoretical mod- 
els declines with increasing halo mass. However, the models 
do not agree in detail at a specific halo mass. And we did 



not find clear evidence that the optical fossil fraction has 
environmental effects in all the theoretical galaxy formation 
models. 

• After applying the bright X-ray luminosity limit in the 
selection of fossil groups, the total fraction of fossil groups 
decreases significantly in small halos. As a consequence, the 
halo mass dependence of the fossil groups seen in optical 
vanishes. 

• From the Lx distribution of groups identified in the 
GIMIC simulations, as shown in Figure [T] we have found ev- 
idence for a clear environmental effect: halos found in lower- 
density regions have relatively lower X-ray luminosity. Their 
median X-ray luminosity can vary by two orders of magni- 
fude, from ~ ]^q38.5 /j^2gj.gg-i f^_2Q- overdensity region in 
GIMIC) to ~ lO*"'^ /i~^ergs~^ (2cr overdensity region in 
GIMIC). 

• When we investigated the properties of central galaxies 
in fossil and non- fossil groups, we have found no differences 
in the magnitude. X-ray luminosity, age, metallicity, concen- 
tration, color or mass-to-light ratio. In particular, the mean 
ages of the central galaxies in fossil groups are not really 
'fossil' in the sense of the word. This casts doubts about 
whether fossil groups are a useful concept to identify a par- 
ticular class of early type galaxies. 

• In addition, we have checked that the satellite galaxies 
number distribution of fossil and non-fossil groups. There is 
not significantly different difference between these two pop- 
ulations as well. 

Our results are cons i stent with the findings of 
Ivon Benda-Beck mann e t al.l (|2008l ). who studied a concor- 
dance ACDM cosmological simulation and pointed out that 
many groups will go through an 'optical fossil phase' that is 
typically ended by renewed infall from the environment. In 
this picture, fossil groups are simply groups that temporarily 
are in a 'fossil phase', such that a significant difference in the 
central galaxy properties can not be expected. It is consis- 
tent with the analysis carried out bv lDariush et all l|201Cf ). 
using also Millennium data, who found that about 90% fos- 
sil groups will become non-fossils after ~ 4Gyr. The other 
important result of our analysis is the absence of strong en- 
vironmental effects at the group scale, at least at the level 
that could be probed with the limited sample sizes that we 
had available in the GIMIC regions. It is well possible that 
there are still weak trends with large-scale overdensity, but 
uncovering those reliably will require larger simulation vol- 
umes. 
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